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Abstract 

We have employed the relativistic coupled cluster theory to calculate the magnetic dipole and 
electric quadrupole hyperfine constants for the ground and low lying excited states of singly ionized 
magnesium. Comparison with experimental and the other theoretical results are done and predictions 
are also made for a few low lying excited states which could be of interest. We have made comparative 
studies of the important many body effects contributing to the hyperfine constants for the different 
states of the ion. 

PACS number(s). : 31.lS.Ar, 31.15.Dv, 31.25.Jf, 32.10.Fn 

1 Introduction 

The hyperfine interactions in alkali metal atoms and alkaline earth ions have been of interest for quite 
a long time pQGj- A number of theoretical investigations including some based on relativistic many- 
body theory have been performed [BlElQi] and they compare reasonably well with experiments. Some 
of the theoretically predicted values could be of experimental interest with the advent of high precision 
techniques involving trapped and laser cooled atoms [J and ions [5] . 

The high precision calculations of different properties of many-electron atoms requires accurate wave- 
functions in the nuclear region as well as the region far from the nucleus. The study of properties like 
hyperfine constants requires the former. Since hyperfine interactions are sensitive to electron correlations, 
the determination of atomic hyperfine constant provides an important test for ab initio atomic structure 
theory 

In this paper we have carried out ab-initio calculations of the magnetic dipole and electric quadrupole 
hyperfine constants and compared the many-body effects for the ground as well as some excited states. 
Section provides the necessary theoretical background to the magnetic dipole (A) and the electric 
quadrupole (B) hyperfine constants. Section gives an overview of the coupled-cluster theory and its 
application to this specific problem. Finally in sectionQ]the results of our calculations are presented and 
discussed. 

2 Theory 

The interaction between the various moments of the nucleus and the electrons of an atom are collectively 
referred to as hyperfine interactions jH]. In this paper we shall consider the interactions between the 
atomic electrons with the nuclear magnetic dipole and electric quadrupole moments. Nuclear spin gives 
rise to a nuclear magnetic dipole moment and the departure from spherical charge distribution in the 
nucleus produces an electric quadrupole moment. 
The hyperfine interaction is given by [2] 

%, = ^MW -T« (1) 
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where and are spherical tensors of rank k, which corresponds to nuclear and electronic parts 

of the interaction respectively. The lowest k — order represents the interaction of the electron with 
the spherical part of the nuclear charge distribution. The eigenstates of the atomic Hamiltonian in the 
presence of a hyperfine interaction are denoted by \IJFMp). Here I and J are the total angular angular 
momentum for the nucleus and the electron state, respectively, and F = I + J with the projection Mp. 

In the first order perturbation theory, the energy corresponding to the hyperfine interaction of the 
fine structure state \ JMj) are the expectation values of H^fs such that 



W[J) = (IJFM F \ J2 M{k) ■ T(k) \UFMf) 

k 

= E(- 1 ) /+J+F { J i I )(I\\M^\\I){J\\T^\\J) 



(2) 



The k = 1 term describes the magnetic dipole coupling of the nuclear magnetic moment with the 
magnetic field created by the electron at the position of the nucleus. The nuclear dipole moment fii is 
defined (in units of Bohr magneton /Uy) as 

= \II) = ( J 7 I \ ) (7|| M« ||7) (3) 

and the operator Tg is given bv[Tfl| 

n i] - E 4 13 - E --Vr r ^ ' (4) 
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Here ~ot is the Dirac matrix and Y£ is the vector spherical harmonics. In Eq.<@J the index j refers to the 
j-th electron of the atom and e is the magnitude of the electronic charge. The magnetic dipole hyperfine 
constant A is defined as 

A= m m , (J||r(1)||J) , ( 5) 

M V I J v /J(J+l)(2J+l)' 1 ; 

and the corresponding magnetic dipole hyperfine energy Wmi is given by 

Wmi=A(I-J) = Aj, (6) 

where K = F(F + 1) - 7(7 + 1) - J {J + 1). 

The second order term in the hyperfine interaction is the electric quadrupole part. The electric 
quadrupole hyperfine constant is defined by putting k — 2 in Eq. IpJ. The nuclear quadrupole moment 
is defined as 

r f-Ef ) = E--7 3 ^ 2) ^)' (7) 



Here, Cq^ — \J (pk+i) ^fcgi with Y^ q being the spherical harmonic. Hence the electric quadrupole hyperfine 
constant B is 



B = 2eQ 



2J(2J-1) 1 J 



(J\\T^\\J), (8) 



(2J+l)(2J + 2)(2J + 3) 
and the corresponding electric quadrupole hyperfine energy Wei is given by 

7£ 3if(if + l)-47(7+l)J(J + l) 
E2 2 27(27- 1)27(27-1) ' K> 

In Eq.(Q]and[3) tq are the single particle reduced matrix element for the electronic part. The reductions 
of the single particle matrix elements into angular factors and radial integral are straightforward by means 
of using the Wigner Eckart theorem. These single particle reduced matrix elements are given by 

(.11 #> |K) = - («|| c« ||*'> (k + k') J dr ( p «Q« + Q* P «') m 

and 
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( K || tf \W) = ( K \\ cf ||«'> J dr ( p ^' + Q^'\ 

where (k\\ Cq \\k') is the reduced matrix element of the spherical tensor and is equal to 
(_l)i+VV(2i + l)(2i' + l)( | * -fi Wfc,0> 



with 



7T(l,k,l')= | 



1 if / + fc + Z' even 
otherwise 



Here the single particle orbitals are expressed in terms of the Dirac spinors with Pj and Qi as large 
and small components respectively. 

3 Overview of the coupled cluster theory : method of calculation 

We start with an N electron closed shell Dirac- Fock (DF) reference state |$). The corresponding corre- 
lated closed shell state is then 

|*) = ex P (T) |$) , (12) 
where T is the core electron excitation operator. Then the Dirac-Coulomb eigenvalue equation is 

ffcxp(T) |$) = Eexp(T) |$) , (13) 
with the Dirac-Coulomb Hamiltonian 



h=j2 («*< ■ + ifi% - i ) mc2 + v N ) + J2 — ■ ( 14 ) 



. ■ Til 

This leads to the exact ground state energy E of the closed-shell part of the system. Here cti and Pi are 
Dirac matrices and Vjv is the nuclear potential. If we consider the DF state |$) as the Fermi vacuum, 
then the normal ordered Hamiltonian is 

H N = H - {<f>\H\<P} = H - E DF . (15) 

If we project ($|exp(— T) from the left we obtain the correlation energy (AE) and if we project 
any of the excited determinant ($*| exp(— T) we additionally get a set of equations which are used to 
obtain the T amplitudes. Using the normal ordered dressed Hamiltonian Hn = exp(— T)H^ exp(T) the 
corresponding equations for correlation energy and amplitudes become 

($|S^|#)=A£7, (16) 

and 

= 0. (17) 

Here the state |<&*) may be singly excited |$^} or double excited \$ r a l) and so on. The indices a,b, - ■ ■ 
refer to holes and p, q, ■ ■ ■ to particles. We have considered the coupled cluster single and double (CCSD) 
approximation, where the cluster operator T is composed of one- and two-body excitation operators, i.e. 
T = 2~i + T2 , and are expressed in second quantization form 

T = T 1+ T 2 =J2 °l a " tP a + \J2 a\a\a h a a e a l. (18) 

ap abpq 

Contracting the ladder operators pi] and rearranging the indices, the amplitude equations can be 
expressed in the form 

A + B{T) -T = 0, (19) 
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where A is a constant vector consisting of the matrix elements ($*| Hn |$), T is the vector of the 
excitation amplitudes and B(T) is the matrix which depends on the cluster amplitudes itself so that Eq. 

i|19H is solved self-consistently. For example, a typical contribution to the term HNT2T2 is 

Kl - ~ E *W£*ft ( 2 °) 

dgrs 

Here Vd grs is the two-electron Coulomb integral and t p ^ d is the cluster amplitude corresponding to a 
simultaneous excitation of two electrons from orbital a and d to p and r respectively. To obtain a full 
set of terms which contribute to this specific excitation, diagrammatic techniques are used. 

The ground state of 25 Mg + contains only one valance electron in the outer most orbital (3s 1/2). To 
calculate the ground state energy of the system we first compute the correlations for the closed shell 
system ( 25 Mg +2 ) using the closed shell coupled cluster approach and then use the technique of electron 
attachment (open shell coupled cluster (OSCC)) method. The energy of the excited state are obtained 
by the same way. In order to add an electron to the kth virtual orbital of the DF reference state we 
define 

K +1 > = 4i$) (21) 

with the particle creation operator aj,. Then by using the excitation operators for both the core and 
valance electron the exact state is defined as [8J: 

K+VexptTHexp^)}!^ 1 ). (22) 

Here {exp(Sfc)} is the normal ordered exponential representing the valance part of the wave operator. 
Here 

Sk = Sik + S 2 k = ^2 a l a kS P k + ^ E a l a \ a b a kS p k q b , (23) 

k^p bpq 

where k stands for valance orbital. Sk contain the particle annihilation operator a^, and because of 
the normal ordering it cannot be connected to any other valance electron excitation operator and then 
{exp(Sfc)} automatically reduces to {1 + Sk}. 
Then we can write the Ea. l22jl as 

K + Vexp(T){l + S fc }|<+ 1 ). (24) 
Following the same procedure as in the closed-shell approach, we obtain a set of equations 

{^ +1 \H^{l + S k }\^ +1 )=H eff (25) 



and 



H N {l + S k }\^ +1 )=H eff U k N+1 {1 + S k }\^ +1 ), (26) 



where the desired roots can be obtained by diagonalizing H e f f . The Eq. l|2fil) is non- linear in Sk because 
the energy difference H e f f is itself a function of Sk ■ Hence, these equations have to solved self-consistently 
to determine the Sk amplitudes. 

Triple excitations are included in our open shell CC amplitude calculations by an approximation that 
is similar in spirit to CCSD(T) [l*2] . The approximate triple excitation amplitude is given by 

opqr VT 2 + ^^2 ( 7] 

°abk — ^ ' ' 

S a + £b + £k — £p — £q — £r 

where iS?'£ are the amplitudes corresponding to the simultaneous excitation of orbitals a, b, k to p, q, r 
respectively and VT and VS are the correlated composites involving V and T, and V and S respectively. 
Efe is the orbital energy of the kth orbital. The above amplitudes (some representative diagrams are given 
in figureOJ are added appropriately in the singles and doubles open shell cluster amplitude equations and 
these equations are then solved self-consistently. We therefore obtain solutions of Siand S2 amplitudes 
taking into consideration the effect of the triple excitations in an approximate way 

The expectation value of any operator O can be written as the normalized form with respect to the 
exact state |v]/ JV + 1 ^ as 
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(y N + 1 \0\V N + 1 ) _ {1 + gt} exp(Tt)Oexp(T) {1 + ff} |$ jv+1 ) 

\ ' - (*JV+i| ^tv+i) - {i + st} exp(Tt) exp(T) {1 + S} l*^ 1 } ' (28) 

For computational simplicity we store only the one-body matrix element of O = exp(TT)0 exp(T). O 
may be expressed in terms of uncontracted single-particle lines The fully contracted part of O will 
not contribute as it cannot be linked with the remaining part of the numerator of the above equation. 

In the LCCSD approximation Eq. l|22l) turns out to be 



\n +1 ) = {l + T + S k }\^ +l ), (29) 

axidH^=H N + H^T. 

The closed and open shell cluster amplitude equations reduce to 

{® k \H N + H^T\<I> } = 0, (30) 

and 



^{l + 5 fc }|< +1 > = K JV+1 S V \^){^\H N \^). (31) 



K N+1 

The orbitals used in the present work are expanded in terms of a finite basis set comprising of Gaussian 
type orbitals (GTO) Q3] 

Fi, k (r) =r fe exp(-a l r 2 ), (32) 

with fc = 0, 1, 2 • • • for s,p, d, ■ ■ ■ type functions, respectively. The exponents are determined by the even 
tempering condition flj] 

a i =ctoP i - 1 . (33) 

The staring point of the computation is the generation of the Dirac-Fock (DF) orbitals jTJj which are 
defined on a radial grid of the form 

n = r [exp(z - l)h - 1] (34) 

with the freedom of choosing the parameters r and h. All DF orbitals are generated using a two 
parameter Fermi nuclear distribution 

9 = 1 + cxp((r - c)/a) ' (35) 

where the parameter c is the half charge radius and a is related to skin thickness, defined as the interval 
of the nuclear thickness in which the nuclear charge density falls from near one to near zero. 

Although we have used a large basis for the generation of the single particle orbitals, the high-lying 
virtual orbitals (above a certain threshold) are kept frozen as their contributions to the high-lying virtuals 
in the T and S amplitudes in the CC equations are negligible. Another advantage of this approximation 
is that it reduces the memory required to store the matrix elements of the dressed operator H and the 
two-electron Coulomb integrals in the main memory, thereby reducing the computational cost. In our 
calculations, we have included all possible single, double and partial triple excitations from the core. 

We have used different ao and (3 for different symmetries. The number of basis functions used to 
generate the even tempered DF states are listed in tableland the values of the parameters ao and (3 
used are also listed. For coupled cluster calculations, we have restricted the basis by imposing an upper 
bound in energy for single particle orbitals and the convergence of our results are shown in table 



4 Results and discussions 

The values of the magnetic dipole hyperfine constant A and electric quadrupole hyperfine constant B 
for different states are given in table El Our calculated value of A for the ground state 3si/ 2 is in good 
agreement (less than 0.6%) with experiment and it is more accurate than a previous calculation (~1%) 
based on second order relativistic many-body perturbation theory (RMBPT) [Hj. This is because unlike 
the previous work our calculation is based on an approach which is equivalent to all-order MBPT. In 
particular, we have taken into account all single, double and a subset of triple excitations to all orders in 
the residual Coulomb interaction. It is, therefore, not surprising that the result of our calculation when 
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S 2 



(a) (b) 

Figure 1: Some typical important diagrams which arise due to the inclusion of triples through En. l(27jl 




Figure 2: Some typical important MBPT diagrams for pair correlation and core-polarization effects. 
The superscripts refer to the order of perturbation and the dashed lines correspond to the Coulomb 
interaction. Particles and holes (labeled by a) are denoted by the lines directed upward and downward 
respectively. The double line represents the O (the hyperfme interaction operator) vertices. The valance 
(labeled by v) and virtual orbitals (labeled by p,q,r..) are depicted by double arrow and single arrow 
respectively, whereas the orbitals denoted by © can either be valance or virtual. 
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Figure 3: The corresponding CC diagrams for pair correlation and core-polarization effects. Here a 
denotes a hole whereas v denotes valance orbital and p,q,r... denote virtual orbitals (particles). 
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Table 1: No. of basis functions used to generate the even tempered Dirac-Fock orbitals and the corre- 
sponding value of ao and /? used. 





Sl/2 


Pl/2 


P3/2 


d 3 /2 


d-5/2 


h/2 


h/2 


Number of basis 


35 


32 


32 


25 


25 


25 


25 


on 


0.00625 


0.00638 


0.00638 


0.00654 


0.00654 


0.00667 


0.00667 


a 


2.03 


2.07 


2.07 


2.19 


2.19 


2.27 


2.27 



Table 2: Value of magnetic dipole (A) and electric quadrupole (B) hyperflne constants in MHz for 
25 Mg + . (T) stands for CC with perturbative partial triples. 



States 


CCSD(T)W 




MBPT 


LCCSD(T) 


(L)CCSD(T) 


Others 1 " 1 ) 


Experiment^ 




A 


B 


A 


A (o) 




A 


A 


3S1/2 


592.86 




602(8) (fe) 


590.73 


597.45 


597.6 


596.2544(5) 






602.46^ 










4si/ 2 


162.32 




164.65^ 


161.79 


163.34 


163.4 




3P1/2 


101.70 




103.20 (c) 


100.69 


102.40 


103.4 




4pi/ 2 


33.83 




34.26 (c) 


33.51 


34.05 






3^3/2 


18.89 


22.91 


19.94^ 


18.95 


19.02 


19.29 




4^3/2 


6.21 


7.48 


6.54( c ) 


6.24 


6.25 






3c! 3 /2 


1.17 


1.26 


1.16^ 


1.16 


1.17 


1.140 




4d 3 / 2 


0.51 


0.48 


0.505 (c) 


0.51 


0.51 







Refs. : (a) : Present work, (6) 0, (c) Present work, (d) : 4 , (e) : j2] 

carried out by using second order RMBPT is in agreement with the result of the earlier calculation (see 
table 0. 

In table[2the results calculated by Safronova et al @] using linearised coupled cluster in singles doubles 
approximation (LCCSD) are given. We have performed LCCSD(T) calculations and the results are given 
in the same table. The reasons for the discrepancy between the two linearised coupled-cluster calculations 
seem to be the inclusion of partial triple excitations by us and the different choice of basis sets. From 
our present calculations it is possible to determine the contributions from the non-linear clusters. These 
contributions vary from 0.3% (3si/ 2 ) to 0.8% (3d 3 / 2 ) for the different states. We find that if we take 
the linear contributions of the T amplitudes but the nonlinear contributions of the combined T and S 
amplitudes and then perform the calculation of the hyperfme constants (EqEHJ a t the linear level (which 
we have named as (L)CCSD(T)), the result (A(3si/ 2 ) — 597.45 MHz) is in excellent agreement with the 
experiment (error is ~0.2%) for the ground state. But theoretically this approach is not complete and 
hence a proper inclusion of the nonlinear effects is desirable. Table 21 shows the contribution of triples 
in the final property. Some partial triples effect are taken into account iteratively according to Eq. I|27JI 
and the S amplitude thus contain some partial triples. The effect of triples in calculation of properties 
(Eq. 12811 thus comes from the triples effect included in S amplitude. Figure shows some of the typical 
diagrams which give rise to the effect due to triples, where figure^ andQ]} correspond to the terms VS 2 
and VT 2 respectively. 

For the present calculations, the number of basis functions actually used are the following : 16si/ 2 , 
14p (p!/ 2 and P3/2), Hd (d 3 / 2 and d 5 / 2 ) and 10/ (/5/ 2 and /V/ 2 ). Excitations from all the core orbitals 
have been considered. 

The important contributions to the magnetic dipole hyperfme constants for different states are given 
in tabled In particular, we have analyzed the contributions from various many body effects and have 
demonstrated that the most important contributions come from core polarization and pair correlation 
effects. The largest contribution comes from O. The next two largest contributions come from (OSi + 
S\0) and(OS' 2 + S 2 0) which correspond to the pair-correlation (PC) and core-polarization (CP) effects 
respectively. The contribution from the corresponding MBPT terms are listed in table EH Figures and 
El represents the pair-correlation and core-polarization diagrams in MBPT and CC respectively. 



8 



Table 3: Convergence of results for A (3si/2 state) using different basis sets 





Sl/2 


Pl/2 


P3/2 


d-3/2 


d-5/2 


/5/2 


Si 12 


A(in MHz) 


Number 


11 


9 


9 


8 


8 


7 


7 


562.10 


of basis 


14 


12 


12 


10 


10 


8 


8 


589.18 




16 


14 


14 


11 


11 


10 


10 


592.86 




17 


15 


15 


11 


11 


10 


10 


592.86 



Table 4: Comparative tables for the value of magnetic dipole (A) hyperfine constants in MHz for 25 Mg + 
calculated using CCSD and CCSD(T). The difference between the values demonstrate the effect of triples. 



States 


CCSD(T) 


CCSD 


Experiment 


3S1/2 


592.86 


593.01 


596.2544(5) 


4si /2 


162.32 


162.39 





Table 5: Comparative study of the contribution from different terms (CCSD(T)) containing the dressed 
operator O in determining the value of magnetic dipole hyperfine constant A for 25 Mg + ior the different 
states. 



Terms 


3si/ 2 


4S1/2 


3^1/2 


4pi/ 2 


3^3/2 


4^3/2 


3^3/2 


4c?3/2 


O 


468.819 


131.616 


77.975 


26.400 


15.337 


5.196 


1.262 


0.563 


OSi + s\o 


40.046 


8.256 


7.344 


2.096 


1.442 


0.412 


0.070 


0.031 


os 2 + s\o 


77.002 


20.587 


14.891 


4.872 


1.832 


0.519 


-0.175 


-0.086 




0.855 


0.129 


0.179 


0.044 


0.035 


0.009 


0.002 


0.0007 


SlOSi 


1.175 


0.162 


0.238 


0.048 


0.0004 


-0.007 


-0.003 


-0.001 


slos 2 
slos 2 


1.175 


0.162 


0.238 


0.048 


0.0004 


-0.007 


-0.003 


-0.001 


5.447 


1.613 


0.997 


0.350 


0.268 


0.101 


0.015 


0.008 



Table 6: Comparative study of the most contributing terms containing the operator O from (CCSD(T)) 
in determining the value of magnetic dipole hyperfine constant A for 25 Mg + for the different states. 



Terms 


3S1/2 


4S1/2 


3pi/ 2 


4pi/ 2 


3^3/2 


4^3/2 


3c?3/2 


4d 3 /2 


O 


468.819 


130.616 


77.975 


26.400 


15.337 


5.196 


1.262 


0.563 


OS! + s\o 


39.713 


8.233 


7.293 


2.093 


1.440 


0.414 


0.070 


0.031 


os 2 + s\o 


77.767 


20.836 


15.153 


4.973 


1.984 


0.569 


-0.194 


-0.096 
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Table 7: Individual contribution from the OSiand OS 2 diagrams for 3si/ 2 state. The values given 
correspond to the respective terms in MHz. 



Orbital OSi OSx OS 2 OS 2 
~3s7/2 2O023 19^357 

lsi /2 13.309 13.196 

2s 1/2 26.665 25.699 

2 Pl/2 0.882 0.867 

2p 3/2 -0.958 -0.878 



We have listed the contributions from the different terms containing the dressed operator O in table 
El and table El gives the contributions from the terms containing the operator O directly. The results 
given in tables El and El show that the CP contribution is larger than the PC in magnitude for all the 
states, although the ratio of the two effects is not uniform. It is important to note that the former 
contribution includes the hyperflne interaction of all the core orbitals while only a specific valence orbital 
is involved in this interaction for the latter (see Figure EJ|. These individual contributions are presented 
in table0for the ground state of Mg + . However, the ground state hyperfine constant A for Ba + exhibits 
exactly the opposite behaviour 221 • Even though Ba + has more core electrons than Mg + , the relativistic 
enhancement of the valence (6s) magnetic dipole hyperfine interaction results in the value of PC exceeding 
that of CP. 



5 Conclusion 

In this paper we have carried out ab-initio relativistic coupled cluster calculations of magnetic dipole (A) 
and electric quadrupole (B) hyperfine constants for the ground and some excited states of 25 Mg + . Some 
contributions from partial triples are also taken into account in our calculation. We have shown that the 
dominant many-body contributions to these properties come from core-polarization and pair-correlation 
effects. 

In addition to comparing with the available experimental data we have also predicted the values of A 
and B for a few states which could be of interest in the future. Using ion-trapping and other experimental 
techniques, it may be possible to measure both the magnetic dipole and electric quadrupole hyperfine 
constants for different states of 25 Mg + , thereby checking the accuracy of our calculations. This would 
constitute an useful test of the validity of the coupled-cluster theory in capturing the many-body effects 
in hyperfine interactions in light atomic systems with a single valence electron. 
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The computation were carried out on our group's 4 CPU E450 Sun Ultra SPARC machine 
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